Feature-based audio content identification

ABSTRACT

An audio signal is sampled and a frequency transform is performed on a succession of sets of samples of the signal to obtain a time dependent power spectrum for the audio signal. Frequency components output by the frequency transform are collected in frequency bands. More than one running average is taken of each semitone frequency band. When the values of two running averages of the same semitone frequency band cross, time information is recorded. Information about average crossing events that have occurred at different times in a set of adjacent semitone frequency bands is combined to form a key. A set of keys obtained from a song provides a means for identifying the song and is stored in a database for use in identifying songs.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application is based upon and claims priority from prior U.S. Provisional Application No. 60/245,799, filed Nov. 3, 2000, the entire disclosure of which is herein incorporated by reference.

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates to signal recognition, and more specifically to a method for automatically identifying audio content such as a sound recording.

2. Description of Related Art

The development of efficient digital encoding methods for audio (e.g., the Motion Picture Experts Group Layer 3 standard known also as MP3), in combination with the advent of the Internet, has opened up the possibility for the entirely electronic sale and distribution of recorded music. This is a potential boon to the recording industry. On the downside, the technical advances also abet the illegal distribution of music. This poses a threat to the propriety interests of recording artists and music distributors. The ease of distributing high fidelity digital copies that do not degrade over successive generations is a far greater problem to the music industry than the limited copying of music onto audio cassettes that occurred prior to the advent of digital audio. Presently, there are a myriad of Internet sites from which a person can obtained bootleg copies of copyrighted music. Thus, for music copyright enforcement, there is a need for a system and method for the automated identification of audio content.

The identification of music from a digital audio file, such as an MP3 file, is not a trivial problem. Different encoding schemes will yield a different bit stream for the same song. Even if the same encoding scheme is used to encode the same song (i.e., sound recording) and create two digital audio files, the files will not necessarily match at the bit level. Various effects can lead to differentiation of the bit stream even though the resulting sound differences as judged by human perception are negligible. These effects include: subtle differences in the overall frequency response of the recording system, digital to analog conversion effects, acoustic environmental effects such as reverb, and slight differences in the recording start time. Further, the bit stream that results from the application of a given encoding scheme will vary depending on the type of audio source. For example, an MP3 file of a song created by encoding the output of a Compact Disc (CD) will not match at the bit level with an MP3 file of the same song created by encoding the output of a stereo receiver.

One solution that has been proposed is to tag copyrighted music by using digital watermarking. Unfortunately numerous methods have been discovered for rendering digital watermarks illegible. In addition, there are forms of noise and distortion that are quite audible to humans, but that do not impede our ability to recognize music. FM broadcasts and audio cassettes both have a lower bandwidth than CD recordings, but are still copied and enjoyed by some listeners. Likewise, many of the MP3 files on the Internet are of relatively low quality, but still proliferate and thus pose a threat to the profitability of the music industry. Furthermore, some intentional evasions of copyright protections schemes involve the intentional alteration or distortion of the music. These distortions include time-stretching and time-compressing. In such cases, not only may the start and stop times be different, but the song durations may be different as well. All such differences may be barely noticeable to humans, but can foil many conventional copyright protection schemes.

Another problem for the music industry and songwriters is the unauthorized use of samples. Samples are short sections of a song that have been clipped and placed into another song. Unless such a sample can be found and identified, the owner of the copyright on the original recording will not be fairly compensated for its use in the derivative work.

There is a need for a method that can identify audio content such as sound recordings despite subtle differences and alterations that arise during processes such as recording, broadcasting, encoding, decoding, transmission, and intentional alteration.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a schematic of an exemplary network on which a preferred embodiment of the present invention can be implemented.

FIG. 2 is a block diagram of software for identifying music according to a preferred embodiment of the present invention.

FIG. 3 is a flow diagram of a process for locating and identifying audio files on a network using an embodiment of the present invention.

FIGS. 4A and 4B are a flow diagram of a process for identifying events from an audio segment according to a preferred embodiment of the present invention.

FIG. 5 is a flow diagram of a process for generating keys from the events produced by the process shown in FIGS. 4A and 4B.

FIG. 6 is a flow diagram of a process for generating keys from the content of a key generator buffer according to a preferred embodiment of the present invention.

FIG. 7 is a flow diagram of a process for filtering percussion events according to a preferred embodiment of the present invention.

FIG. 8 is a flow diagram of a process for using keys to compare two audio segments.

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS

Objects, features, and advantages of the present invention will become apparent from the following detailed description. It should be understood, however, that the detailed description and specific examples, while indicating preferred embodiments of the present invention, are given by way of illustration only and various modifications may naturally be performed without deviating from the present invention.

FIG. 1 is a schematic of an exemplary network 100 on which the present invention can be implemented. The network includes a first web server 102 and a second web server 104 communicatively coupled to a network 106 such as the Internet through a first bi-directional data link 108 and a second bi-directional link 110, respectively. The first and second web servers 102 and 104 have file storage memories such as hard drives for storing files including audio files that can be downloaded by clients.

A first client computer 112 and a second client computer 114 are communicatively coupled to the network 106 through a third bi-directional data link 116 and fourth bi-directional data link 118, respectively. The first and second client computers are capable of downloading files including audio files from the first and second web servers. The network components described thus far describe a system that can be used for the distribution of music. Audio files storing sound recordings of music are stored on the web servers 102 and 104.

A copyright enforcement client computer (CECC) 120 is communicatively coupled to the network 106 through a fifth bi-directional data link 122. A computer readable memory medium 124 such as a CD-ROM is provided for loading software onto the CECC for carrying out methods such as those described in detail below. For example, the software can operate to spider over (i.e., search) web sites such as web sites hosted by the first and second web servers 102 and 104, retrieve audio files included on the web sites, and identify the audio files (e.g., identifies songs included in the audio files).

The first and second web servers 102 and 104, the client computers 112 and 114, and the CECC can be any conventional computer systems such as IBM PC-compatible computers. As is known, each IBM PC-compatible computer can include a microprocessor, basic input/output system read-only memory (BIOS ROM), random access memory (RAM), hard disk drive storage, removable computer readable medium storage (e.g., a CD-ROM drive), a video display adapter card, a video monitor, a network interface (e.g., modem), a keyboard, a pointing device (e.g., mouse), a sound card, and speakers.

The first through fifth bi-directional data links 108, 110, 116, 118, and 122 may include Digital Subscriber Lines (DSL), T1 lines, or dial-up modem connections. The first and second web servers 102 and 104 can be provided with Apache Server web server software (produced by Apache Software Foundation of Lincoln, Nebr.) running under a UNIX operating system. The first and second client computers 114 and 116 can be loaded with web browsers such as Netscape Navigator (produced by America Online of Dulles, Va.) running under a Windows operating system (produced by Microsoft Corporation of Redmond, Wash.). The web browser preferably operates along with a plug-in application for decoding audio files and providing an audio signal to the sound card. Alternatively, a separate application or dedicated hardware can be provided for decoding audio files for playback. The two web servers 102 and 104, the two client computers 112 and 114, and the CECC are also loaded with communication protocol stack software so as to be able to establish network connections such as TCP/IP connections. Further, the CECC is loaded with an operating system (such as Windows or UNIX) and a CECC application, such as one that operates as described below with reference to FIGS. 2-8.

FIG. 2 is a block diagram 200 of software for identifying audio content according to a preferred embodiment of the present invention. A web server application 202 includes a file system 204 that includes one or more audio files (e.g., in MP3 format). The web server application 202 operates along with a server side communication protocol stack 206 that is capable of supporting network (e.g., TCP/IP) connections. The web server application 202 and server side communication protocol stack 206 are run on the first and second web servers 102 and 104. A TCP/IP connection 208 is established between the server communication protocol stack 206 and a client communication protocol stack 210. The client communication protocols stack is utilized by a copyright enforcement application 212.

The copyright enforcement application includes a web spider module 214 and a content identifier module 216. The web spider module 214 searches web sites for audio files. When an audio file is located, it is downloaded for analysis. The content identifier module 216 receives audio files from the web spider module 214 and decodes them to obtain audio signal information. The content identifier module includes a key database that contains numerous keys derived from numerous songs. For each song (i.e., sound recording unit) in the database, there is a set of keys. The set of keys provides a means of identifying a song or a segment of a song. A section of a song will have a corresponding subset of keys that allow the section to be identified, thus retaining the ability to identify the presence of only a portion of a song.

The key database preferably takes the form of a key table. Each row of the key table includes a key sequence in a first column, a time group value in a second column, and a song ID (e.g., title) in a third column. The key database is constructed by applying a key generation program to known songs, and associating each key obtained by the key generation program with the title in the key database. The time group is a time (measured in units of a predetermined interval relative to the start of the song) at which an audio feature from which a key is derived occurred. Each time group includes events (explained further below) that occurred during a period of time equal to the predetermined interval. Two different recordings of the same song, especially in the case of bootlegged copies, can have slightly different start times (e.g, an audio file may start a few seconds before or after the actual beginning of the song). This variance presents a difficulty in comparing two audio files that the present invention surmounts.

FIG. 3 is a flow diagram of a process 300 that is performed by the copyright enforcement module 210 to locate and identify audio files on the Internet (e.g., web or FTP sites) according to one embodiment of the present invention. In step 302, the web spider module 214 spiders a portion of the Internet and finds a number of audio files. One method for identifying audio files is to compare the file extensions of Uniform Resource Identifiers (URI) found within web pages to a list of known audio file extensions. In step 304, the web spider module 214 downloads the audio files found in step 302.

In step 306, the content identifier module 218 analyzes the audio files that were downloaded in step 304 using an algorithm that generates a unique signature (preferably a set of keys) for each audio file. In step 308, the unique signature is used to look up information such as the title of each audio file that is identified in a database of copyrighted sound recordings (e.g., songs). In step 310, an identification of the copyrighted song (e.g., it's title) and the URI at which it was found on the web are recorded and/or transmitted to a third party (e.g., the owner of the copyright). The list of URIs corresponding to a given copyrighted work can thus be compiled for use in copyright enforcement actions.

FIG. 4A is a first part of a flow diagram of a process 400 for generating an event stream for an audio segment (e.g., song) according to a preferred embodiment of the present invention. The process 400 accepts an audio signal as input and outputs sequence of “events”. In some embodiments of the present invention, the audio signal is reproduced from an MP3 file. In step 402, an audio signal is sampled. In one embodiment, the audio signal is sampled at about 22050 Hz or lower. This allows frequency components up to 11 KHz to be accurately determined. It is advantageous to use an upper frequency limit of about 11 KHz because 11 KHz is about the frequency cutoff for FM broadcast radio, and it is desirable to be able to generate the same set of keys for a song regardless of whether the song recording was at one point transmitted through FM radio or obtained directly from a high quality source (e.g., a CD).

In step 404, for each successive test period the spectrum of the audio signal is computed. The duration of the test period preferably ranges from about {fraction (1/43)} of a second to about {fraction (1/10.75)} of a second, and more preferably the test period is about {fraction (1/21.5)} of a second. The spectrum of the audio signal is preferably analyzed using a fast Fourier transform (FFT) algorithm. The accuracy of spectrum information obtained using an FFT algorithm can be improved by averaging together the results obtained by applying the FFT to several successive periods (sample sets). In preferred embodiments of the present invention, spectrum information is improved by averaging together the results obtained by applying the FFT to two or more successive periods, and preferably 3 or more successive periods, and even more preferably 4 successive periods. According to one exemplary embodiment of the present invention, the spectrum associated with a given test period having a duration of {fraction (1/21.5)} of a second is obtained by sampling an audio signal at a rate of 22050 Hz and averaging together the results obtained by applying an FFT algorithm to four successive periods, each of which has a duration of {fraction (2/21.5)} seconds and includes 2048 samples.

Step 404 can be accomplished by using an FFT algorithm run on the microprocessor of the CECC 120. Alternatively, the CECC could be provided with FFT hardware for performing step 404. Other spectrum analyzers, such as a filter bank, can alternatively be used for carrying out step 404. Additionally, in process 404, successive sets of samples can alternatively be projected onto another type of basis besides a Fourier basis. One particular alternative to the Fourier basis is a wavelet basis. Like Fourier basis functions, wavelets are also localized in the frequency domain (although to a lesser degree). Wavelets have the added property that they are localized in the time domain as well. This opens up the possibility of projecting the audio signal as a whole, rather than successive sample sets of the audio signal onto a wavelet basis, and obtaining time dependent frequency information about the signal.

One common set of frequencies used in composing music are the notes of the even-tempered scale. The even tempered scale includes notes that are equally spaced on a logarithmic scale. Each note covers a frequency band called a “semitone”. The inventors have determined that improved signatures can be obtained by collecting spectral power in discrete semitone bands as opposed to the evenly spaced frequency bands output by an FFT algorithm. In step 406, the spectrum information (e.g., Fourier frequency components) obtained in step 404 are collected into a number of semitone frequency bands or channels.

In step 408, a first average of the power in each semitone frequency channel is taken over the last T1 seconds. In step 410, a second average of the power in each semitone frequency channel is taken over the last T2 seconds, where T2 is greater than T1. T1 is preferably from about {fraction (1/10)} to about 1 second. T2 is preferably larger than T1 by a factor of from 2 to 8. According to a one exemplary embodiment of the present invention, T2 is equal to one second, and T1 is equal to one-quarter of a second. The “events” mentioned above occur when the value of the first average crosses the second average.

In step 412, the values of the first and second averages are recorded for each semitone channel. Recording is done so that it can be determined during the following test period whether the first average crossed the second average. In step 414, for each semitone channel it is determined if the first average crossed the second average. This is done by comparing the inequality relation between the first and second averages during the current test period to the inequality relation for the last period. Although comparison between only two averages has been discussed above, it is possible according to alternative embodiments of the present invention to use more than two averages, and identify events as the crossing points between different sub-combinations of the more than two averages.

In the vicinity of an extremum (local maximum or minimum) in a semitone frequency channel, the two averages will cross. Rather than looking for the crossing point of two running averages with different averaging periods, another type of peak detector (e.g., an electronic circuit) could be used. Such could advantageously be used in combination with an FFT in an implementation of the present invention that is implemented predominately in hardware, as opposed to software.

Rather than looking for a peak in the signal in a frequency channel, another type of curve characteristic such as an inflection point could be used as a trigger event. An inflection point can be found by calculating a second derivative of a frequency channel by operating on three successive values of the power in a given frequency channel, and identifying a time at which the second derivative changes from positive to negative or vice versa. The second derivative can be approximated using function (time dependent frequency component) values for three successive points in the following formula.

(F(N+2)−2F(N+1)+F(N))/ΔT ²

where F(I) is the value of the function at the i^(th) time (e.g., at the i^(th) test period), and ΔT is the interval between successive function values (e.g., the duration of the test period).

At an extremum of a time dependent frequency component, its first derivative is equal to zero. At an inflection point of a time dependent frequency component, its second derivative is equal to zero. Extrema and inflection points are both kinds of events. More generally events can be defined as points (i.e., points in time) at which an equation involving a time dependent frequency component derivative of one or more orders of the time dependent frequency components, and/or integrals involving the time dependent frequency components is satisfied. To allow their use in identifying different audio content, an essential part of the definition of “events” is that they occur at a subset of test periods, not at each test period.

Step 416 is a decision block, the outcome of which depends on whether averages for a semitone channel crossed. Step 416 is tested for each semitone channel. If averages for a semitone channel were not crossed during the current test period, then in step 418 it is determined if the audio signal is over. If the audio stream is finished, then the process 400 terminates. If the audio signal is not finished, then the process 400 is advanced to the next test period and the process continues with step 404. If on the other hand, averages did cross during the last test period then the process 400 continues with step 422 in which each event is assigned to the current time group and information related to the average crossing event is generated.

Event information preferably includes the time group for the event, the test period for the event, the semitone frequency band of the event, and the value of the fast average (average over T1) at the time of crossing. Event information can be recorded in a memory or storage device associated with the CECC. Each time group covers a period of time that is longer than a test period, and preferably time groups cover successive periods of time equal to from ¼ to 2 seconds, and more preferably each time group covers a period of from one-half to three-quarters of a second. Grouping events into successive time groups has the advantage that keys obtained by processing two recordings of the same song will tend to match more completely despite the fact that one or both of the recordings may have some distortions (e.g., distortions that arise in the course of recording on magnetic tape).

In step 424, the process 400 is incremented to the next test period. In step 426, it is determined if the audio segment (e.g., song) is finished. If the audio segment is finished then the process 400 terminates. If the audio segment is not finished, then the test period is incremented and the process loops back to step 404.

Thus, the result of the process is to take an audio signal and produce a plurality of events. Each event is assigned to a semitone frequency band in which it occurred and a time group (interval) within which it occurred. The events can be stored in a memory (e.g., RAM in the CECC 120). The events can be stored in a buffer from which they are successively read by one or more key generator processes. The events output by the process could be in the form of an event stream, which is to say that after each time group, all the events occurring within the time group could be written to memory and thereby made available for further processing. An alternative is to write all the events for a song to memory or storage at one time.

FIG. 5 is a flow diagram of a key generator process for generating keys from the events produced by a process such as that shown in FIGS. 4A and 4B. The events output by process 400 are processed by a plurality of key generator processes 500. Each of the plurality of key generator processes is assigned to one semitone frequency band that is designated as its main frequency. However, each key generator also uses events that occur in other semitone frequency bands near its main frequency. Preferably each key generator monitors from 5 to 15 semitone frequency bands. If the number of frequency bands monitored is too few, the resulting keys will not be as strongly characteristic of the particular audio segment. On the other hand, a higher number of frequency bands will result in higher computational expense for computing and comparing keys, greater memory requirements for storing keys, and potential performance loss due to key saturation in the key table from the increased number of keys. According to one embodiment of the present invention, each key generator monitors its main semitone frequency band and four other semitone frequency bands, two on each side of the main semitone frequency band.

Referring now to FIG. 5, in step 502 each successive time group of events output by process 400 is monitored for events occurring within the semitone frequency bands assigned to this key generator. Step 504 is a decision block, the outcome of which depends on whether the key generator detected (e.g., by reading from memory) any new events in step 502. If not, then in step 514, the process 500 is incremented to the next time group and loops back to step 502. If, on the other hand, new events did occur in the time group and semitone frequency bands checked, then in step 506 the new events are written to a key generator buffer for the key generator under consideration, and the events for the oldest time group that were stored in the key generator buffer are deleted. In one exemplary embodiment, the buffer can be seen as an array in which the rows correspond to time groups and the columns to frequency bands. Thus, in the embodiment of the present invention mentioned above, there would be five columns for each of the semitone frequency bands monitored by each key generator.

The key generator buffer preferably includes events from 3 to 7 time groups. More preferably, events from five or six time groups are maintained in each key buffer array. Note that in this embodiment not all time groups are represented in the key generator buffer. As shown in FIG. 5, if no events occur in the semitone frequency bands for a key generator in a certain time group, then no change will be made to the key generator buffer. In other words, a blank row will not be recorded. Therefore, each time group recorded in the key generator buffer includes at least one event.

Step 508 is decision block whose outcome depends on whether an event that occurred in the current time group (e.g., current pass through program loop) is a trigger event. According to a preferred embodiment of the present invention, a trigger event is an event that occurs at the main frequency assigned to this key generator. If a trigger event did not occur, then the process loops back to step 514. If a trigger event did occur, then the process continues with step 510 in which keys are generated from the contents of the key generator buffer. The process 500 continues until all of the events produced by process 400 have been processed.

FIG. 6 is a flow diagram of a process for generating keys from the contents of a key generator buffer according to one embodiment of the present invention. In particular, the process 600 shows in detail one embodiment of the implementation of step 510 of FIG. 5. In step 602, for each key generator (as explained above there are a plurality of key generators carrying out process 500) and for each trigger event for the key generator under consideration, one or more different combinations of events from the key generator buffer are selected. Each combination includes only one event from each time group. (There may be more than one event for each time group in each key generator buffer.) According to a preferred embodiment of the present invention, not all possible combinations are selected, rather only combinations for which a power associated with each event changes monotonically from one event to the next in the combination are selected.

In this embodiment, the order of events within a combination corresponds to the time group order. The power associated with each event is preferably the magnitude of the fast (first) average at the test period at which the event occurred. In this embodiment, less than all of the possible combinations of keys will be taken, so that the total number of keys for a given audio segment will tend to be reduced which leads to lower memory and processing power requirements. On the other hand, there will be enough keys that the identity of the song will be well characterized by (i.e., strongly correlated to) the set of keys generated from the song. According to an alternative embodiment, only a single combination is selected from the contents of the key generator buffer. The single combination includes the event associated with the highest fast average power from each time group. According to another alternative embodiment, all the different combinations of events taking one event from each time group are taken.

In step 604, for each selected combination of events a key sequence is composed that includes a sequence of numerical values of frequency offsets (relative to the main key generator frequency) for the sequence of events from each combination formed in step 602. Each frequency offset is the difference between the frequency of the semitone band in which the event occurred and the main frequency of the key generator. In step 606, test period information (e.g., a sequence number for the test period of the trigger event, where the sequence number for the first test period for each song is designated by the number one) for the trigger event is associated with the key sequence.

In step 608, the key which includes the key sequence and the test period information is associated with a song (or other audio) identifier or ID (e.g., title). Process 600 includes step 608 in the case that known songs are being used to construct a song database against which unknown songs will be compared. In comparing two songs, both the key sequence and test period information will be used, as described further below with reference to FIG. 8. The song database can take the form of a table including three columns and a plurality of rows. The first column includes key sequences, the next column includes corresponding test periods associated with the key sequences, and the final column includes an identification of the song from which the keys in the row were obtained.

While the processes described above can be used to identify audio content, it is advantageous to filter percussion events. More specifically, percussion sounds in a song, if not filtered, typically account for high percentage of the events output by process 400. In the interest of saving computer resources (e.g., memory and processing power) and obtaining a more characteristic set of keys, it is desirable to reduce the number of percussion events such as by eliminating some percussion events before events are processed by the key generator process 500. It has been recognized by the inventors that percussion sounds lead to events being triggered during the same test period in adjacent semitone frequency bands. For example, percussion sounds can lead to events occurring in a sequence of 2 or more adjacent semitone frequency bands.

FIG. 7 is a flow diagram of a process used in a preferred embodiment to filter percussion events from the events produced by the process of FIGS. 4A and 4B. In step 702, for each successive test period it is determined if multiple events occurred in a sequence of two or more adjacent semitone frequency bands. Step 704 is a decision block, the outcome of which depends on whether multiple events in adjacent frequency bands occurred. A threshold of some predetermined number of events occurring in adjacent frequency bands used in the process. Preferably, a lower limit on the number of adjacent frequency bands in which events must be found (in order to consider that the events were produced by a percussion sound) is set at three or more. According to an exemplary embodiment of the present invention, events must occur in three successive semitone frequency bands for the outcome of step 704 to be positive.

If the outcome of step 704 is negative, then the process continues with step 708 in which the process increments to the next test period and loops back to step 702. If, on the other hand, the outcome of step 704 is positive, then the process 700 continues with step 706 in which each sequence of events that occurred during the same test period in adjacent frequency bands is pared down to a single event. All of the events except for the event in the sequence that has the highest fast average value are deleted from the event stream produced by process 400. Alternatively, instead of deleting all but one, up to a certain predetermined number of events can be retained.

The processes described above produce keys for a sound recording based on the features (i.e., events) contained in the sound recording. Thus, the processes can be ran on known audio content to construct a feature database of the known audio content during a storage phase. After the database is created, during a retrieval phase the above processes can be used to extract features from unknown audio content and then the database can be accessed to identify the audio content based on the features that are extracted. For example, the same processes can be run on the unknown audio content to extract features in real time (or even faster), and then the audio content is identified with the best match in the database. In one embodiment, a best match can be reported for each predetermined interval (e.g., 10 to 30 seconds) of the audio content.

FIG. 8 is a flow diagram of a song identification process that uses the keys generated in the processes of FIGS. 5 and 6 to identify an audio segments. A song database (such as that described above) is used to identify an unknown song such as a song downloaded from a web site in step 304 of process 300. The key sequence field (column) of the song database can be used as a database key. The records (rows) of the song database are preferably stored in a hash table for direct lookup. The identification process 800 is an exemplary implementation of step 308 of FIG. 3.

In step 802, keys are generated from a song to be identified (for example, by carrying out the processes shown in FIGS. 5 and 6). In step 804, each key in the set of keys generated in step 804 is looked up in a song database that includes keys for a plurality of songs. The key sequence part (as opposed to the test period part) of each key is used as a database key. In other words, the song database is searched for any entries that have the same key sequence as a key sequence belonging to a key obtained from the song to be identified. More than one key in the song database can have the same key sequence, and furthermore by happenstance more than one song in the song database can share the same key sequence. In step 806, for each key in the database that matched (by key sequence) one or more keys in the song database, an offset is calculated by taking the difference between a test period associated with the key being looked up and a test period associated with each matching key in the song database.

In step 808, the offsets are collected into offset time groups. The offset time groups for the offsets are distinct from the time groups used in key generation. According to a preferred embodiment, an offset time group will be equal to from 2 to 10 test periods. By way of illustration, if each offset time group were 5, then any pair of keys for which the difference determined in step 806 was between 0 and 5 would be assigned to a first offset time group, and any pair of keys for which the difference was between 6 and 10 would be assigned to a second offset time group. According to an exemplary embodiment of the present invention, each offset time group is equal to 5 test periods.

In step 810, for each song that has keys that match keys in the song to be identified, and for each offset time group value that was determined in step 808 and involved keys for a given song in the song database, a count is made of the number of matching keys that had the same time group offset value. One can visualize step 810 in the following way, which may also be used as a basis for an implementation approach. A temporary table is constructed where each row corresponds to a song from the song database that had one or more key matches with the song to be identified. The first column includes names of the songs. In the second column, adjacent each song name there is a value of the offset time group that was found between keys found for the named song in the song database, and matching keys from the song to be identified. After completing step 810, the third column will include counts of the number of key matches corresponding to a particular song identified in the first column, that had the same offset time group as identified in the second column. The table might appear as follows.

TABLE 1 COUNT OF KEY SEQUENCE MATCHES OFFSET VALUE FOR THIS SONG AND (UNITS OF TIME WITH THIS OFFSET SONG TITLE GROUP INTERVAL) VALUE Title1 3 1 Title1 4 1 Title2 2 2 Title2 3 107 Title3 5 1 Title2 8 1

If the song to be identified is in the database, then one particular time group offset value will accumulate a high count. In other words, a high number of matching pairs of keys will be found to have some particular value of offset time group. In the example above, the song entitled Title2 has a count of 107 for an offset time group of 3. For example, the time group offset may arise because the specific recording that is being identified started a few seconds after the recording of the song used to generate keys for the song database, or because a small segment of the song is being identified.

In step 812, the song from the song database that has the highest count of matching keys with the same offset is identified. In decision block 814, the count is compared to a threshold value. The threshold can be set based on the particular application or through a determination of the minimum value for the highest counts that are found when songs actually match, and the maximum value of the highest counts when songs tested do not match any songs in the database. The value of the threshold used also depends on the specific embodiment chosen for step 602 discussed above, as this determines the total number of keys.

Rather than comparing the count to a threshold, it is possible instead in step 812 to compare a threshold with the ratio of the highest count to the total number of keys generated from the song to be identified. Another alternative is to compare a threshold with the ratio of the highest count to the average of the remaining counts. These latter two alternatives can also be viewed as comparing the highest count to a threshold, although in these cases the threshold is not fixed. If, as would be the case when the song to be identified is not in the database, the count does not meet the threshold criteria, then the song identification process 800 terminates. Additional steps may be provided for reporting (e.g., to a user) that the song to be identified could not be identified. If on the other hand the count does meet the threshold criteria, then in step 814 information identifying the song that had the highest count (which met the threshold criteria) is output. In further embodiments, the processes of the present invention are used to identify segments of songs.

Accordingly, the present invention provides methods for identifying songs that are robust in terms of their ability to handle distortion and alteration. Furthermore, the method is also efficient in terms of computational complexity and memory requirements. In preferred embodiments, the processes for generating an event stream, filtering percussion events, generating keys, and looking up the keys in a song database are conducted in real time. The computational expense of these processes in the preferred embodiments is low enough to allow them to run in real time on a common personal computer.

The present invention can be realized in hardware, software, or a combination of hardware and software. Any kind of computer system—or other apparatus adapted for carrying out the methods described herein—is suited. A typical combination of hardware and software could be a general purpose computer system with a computer program that, when being loaded and executed, controls the computer system such that it carries out the methods described herein.

The present invention can also be embedded in a computer program product, which includes all the features enabling the implementation of the methods described herein, and which—when loaded in a computer system—is able to carry out these methods. Computer program means or computer program in the present context mean any expression, in any language, code or notation, of a set of instructions intended to cause a system having an information processing capability to perform a particular function either directly or after either or both of the following a) conversion to another language, code or, notation; and b) reproduction in a different material form.

Each computer system may include, inter alia, one or more computers and at least a computer readable medium allowing a computer to read data, instructions, messages or message packets, and other computer readable information from the computer readable medium. The computer readable medium may include non-volatile memory, such as ROM, Flash memory, disk drive memory, CD-ROM, and other permanent storage. Additionally, a computer medium may include, for example, volatile storage such as RAM, buffers, cache memory, and network circuits. Furthermore, the computer readable medium may include computer readable information in a transitory state medium such as a network link and/or a network interface, including a wired network or a wireless network, that allow a computer to read such computer readable information.

While there has been illustrated and described what are presently considered to be the preferred embodiments of the present invention, it will be understood by those skilled in the art that various other modifications may be made, and equivalents may be substituted, without departing from the true scope of the invention. Additionally, many modifications may be made to adapt a particular situation to the teachings of the present invention without departing from the central inventive concept described herein. Therefore, it is intended that the present invention not be limited to the particular embodiments disclosed, but that the invention include all embodiments falling within the scope of the appended claims. 

What is claimed is:
 1. A method for identifying audio content, said method comprising the steps of: obtaining an audio signal; analyzing the power spectrum of the audio signal so as to obtain a plurality of time dependent frequency components; and detecting a plurality of events, each of the events being a crossing of the value of a first running average and the value of a second running average, wherein the first running average is an average over a first averaging period of a first subset of the time dependent frequency components, and the second running average is an average over a second averaging period, which is different than the first averaging period, of the first subset of the time dependent frequency components.
 2. The method according to claim 1, further comprising the steps of: detecting a set of events occurring at nearby times in adjacent frequencies bands; and forming a key by combining at least a subset of the set of events.
 3. The method according to claim 1, further comprising the step of determining a time dependent frequency component power corresponding to each event.
 4. The method according to claim 1, wherein the analyzing step includes the sub-steps of: sampling the audio signal to obtain a plurality of audio signal samples; taking a plurality of subsets from the plurality of audio signal samples; and performing a Fourier transform on each of the plurality of subsets to obtain a set of Fourier frequency components.
 5. The method according to claim 4, wherein the analyzing step further includes the sub-step of averaging together corresponding Fourier frequency components obtained from two or more successive subsets selected from the plurality of subsets.
 6. The method according to claim 5, wherein the analyzing step further includes the sub-step of collecting Fourier frequency components into a plurality of semitone frequency bands.
 7. The method according to claim 1, wherein the detecting step includes the sub-steps of: keeping the first running average over the first averaging period of a first subset of the time dependent frequency components so as to obtain a first series of averages for the first averaging period; keeping the second running average over the second averaging period of the first subset of the time dependent frequency components so as to obtain a second series of averages for the first averaging period; and recording a plurality of event times, each of the event times being a time at which there occurs one of the detected events of the first running average crossing the second running average.
 8. The method according to claim 1, wherein the first averaging period is between approximately {fraction (1/10)} of a second and approximately 1 second, and the second averaging period is from approximately 2 to 8 times as long as the first averaging period.
 9. The method according to claim 1, further comprising the step of collecting the plurality of events in a plurality of time groups each of which covers an interval of time.
 10. The method according to claim 9, further comprising the step of: in response to detecting each event in each of the plurality of time dependent frequency components, selecting one or more combinations of events from a plurality of events that occurred within a number of time groups, and within a number of time dependent frequency components.
 11. The method according to claim 10, wherein the selecting step includes the sub-step of selecting one or more combinations of events from a plurality of events that occurred within a number of time groups, and within a number of time dependent frequency components, taking only one event at a time from each time group.
 12. The method according to claim 10, further comprising the step of forming a plurality of keys from the one or more combinations each of which comprises a time to be associated with the combination of events, and a key sequence including information about each event in the combination.
 13. A method for forming an identifying feature of a portion of a recording of audio signals, said method comprising the steps of: performing a Fourier transformation of the audio signals of the portion into a time series of audio power dissipated over a first plurality of frequencies; grouping the frequencies into a smaller second plurality of bands that each include a range of neighboring frequencies; detecting power dissipation events in each of the bands; and grouping together the power dissipation events from mutually adjacent bands at a selected moment so as to form the identifying feature, wherein each of the power dissipation events is a crossing of the value of a first running average and the value of a second running average, the first running average is an average over a first averaging period of the audio power dissipated, and the second running average is an average over a second averaging period, which is different than the first averaging period, of the audio power dissipated.
 14. The method according to claim 13, further comprising the step of integrating power dissipation in each of the bands over a predetermined period.
 15. The method according to claim 14, wherein each of the power dissipation events is a crossover of rolling energy dissipation levels over time periods of different lengths.
 16. A method of determining whether an audio stream includes at least a portion of a known recording of audio signals, said method comprising the steps of: forming at least a first identifying feature based on the portion of the known recording and at least a second identifying feature based on a portion of the audio stream using the method of claim 13; storing the first identifying feature in a database; and comparing the first and second identifying features to determine whether there is at least a selected degree of similarity.
 17. The method according to claim 16, wherein each of the power dissipation events is a crossover of rolling energy dissipation levels over time periods of different lengths.
 18. A computer-readable medium encoded with a program for identifying audio content, said program containing instructions for performing the steps of: obtaining an audio signal; analyzing the power spectrum of the audio signal so as to obtain a plurality of time dependent frequency components; and detecting a plurality of events, each of the events being a crossing of the value of a first running average and the value of a second running average, wherein the first running average is an average over a first averaging period of a first subset of the time dependent frequency components, and the second running average is an average over a second averaging period, which is different than the first averaging period, of the first subset of the time dependent frequency components.
 19. The computer-readable medium according to claim 18, wherein said program further contains instructions for performing the steps of: detecting a set of events occurring at nearby times in adjacent frequencies bands; and forming a key by combining at least a subset of the set of events.
 20. The computer-readable medium according to claim 18, wherein the analyzing step includes the sub-steps of: sampling the audio signal to obtain a plurality of audio signal samples; taking a plurality of subsets from the plurality of audio signal samples; and performing a Fourier transform on each of the plurality of subsets to obtain a set of Fourier frequency components.
 21. The computer-readable medium according to claim 18, wherein the detecting step includes the sub-steps of: keeping the first running average over the first averaging period of a first subset of the time dependent frequency components so as to obtain a first series of averages for the first averaging period; keeping the second running average over the second averaging period of the first subset of the time dependent frequency components so as to obtain a second series of averages for the first averaging period; and recording a plurality of event times, each of the event times being a time at which there occurs one of the detected events of the first running average crossing the second running average.
 22. A computer-readable medium encoded with a program for forming an identifying feature of a portion of a recording of audio signals, said program containing instructions for performing the steps of: performing a Fourier transformation of the audio signals of the portion into a time series of audio power dissipated over a first plurality of frequencies; grouping the frequencies into a smaller second plurality of bands that each include a range of neighboring frequencies; detecting power dissipation events in each of the bands; and grouping together the power dissipation events from mutually adjacent bands at a selected moment so as to form the identifying feature, wherein each of the power dissipation events is a crossing of the value of a first running average and the value of a second running average, the first running average is an average over a first averaging period of the audio power dissipated, and the second running average is an average over a second averaging period, which is different than the first averaging period, of the audio power dissipated.
 23. A system for identifying a recording of an audio signal, said system comprising: an interface for receiving an audio signal to be identified; a spectrum analyzer for analyzing the power spectrum of the audio signal so as to produce a plurality of time dependent frequency components from the audio signal; an event detector for detecting a plurality of events in each of the time dependent frequency components; and a key generator for grouping the plurality of events by frequency and time, and assembling a plurality of keys based on the plurality of events, wherein each of the events detected by the event detector is a crossing of the value of a first running average and the value of a second running average, the first running average is an average over a first averaging period of a first subset of the time dependent frequency components, and the second running average is an average over a second averaging period, which is different than the first averaging period, of the first subset of the time dependent frequency components.
 24. The system according to claim 23, wherein the event detector is a peak detector.
 25. The system according to claim 23, further comprising a database of keys of known recordings of audio signals.
 26. A system for forming an identifying feature of a portion of a recording of audio signals, said system comprising: means for performing a Fourier transformation of the audio signals of the portion into a time series of audio power dissipated over a first plurality of frequencies; means for grouping the frequencies into a smaller second plurality of bands that each include a range of neighboring frequencies; means for detecting power dissipation events in each of the bands; and means for grouping together the power dissipation events from mutually adjacent bands at a selected moment so as to form the identifying feature, wherein each of the power dissipation events detected by the means for detecting is a crossing of the value of a first running average and the value of a second running average, the first running average is an average over a first averaging period of the audio power dissipated, and the second running average is an average over a second averaging period, which is different than the first averaging period, of the audio power dissipated. 